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Stochastic background of gravitational waves emitted by 
magnetars 
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ABSTRACT 

Two classes of high energy sources in our galaxy are believed to host magnetars, 
neutron stars whose emission results from the dissipation of their magnetic field. The 
extremely high magnetic field of magnetars distorts their shape, and causes the emis- 
sion of a conspicuous gravitational waves signal if rotation is fast and takes place 
around a different axis than the symmetry axis of the magnetic distortion. Based on 
a numerical model of the cosmic star formation history, we derive the cosmological 
background of gravitational waves produced by magnetars, when they are very young 
and fast spinning. We adopt different models for the configuration and strength of the 
internal magnetic field (which determines the distortion) as well as different values of 
the external dipolc field strength (which governs the spin evolution of magnetars over 
a wide range of parameters). We find that the expected gravitational wave background 
differs considerably from one model to another. The strongest signals are generated 
for magnetars with very intense toroidal internal fields (~ 10 16 G range) and external 
dipole fields of ~ 10 14 , as envisaged in models aimed at explaining the properties of 
the Dec 2004 giant flare from SGR 1806-20. Such signals should be easily detectable 
with third generation ground based interferometers such as the Einstein Telescope. 
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1 INTRODUCTION 

It is well known that a variety of astrophysical processes are 
able to generate a stochastic gravitational wave background 
( GWB), with distinct spectral propert i es and features 
Ferrari et all 1 1999al; fFerrari et all Il999bl; [Schneider et all 



l200d . l200ll : lRegimbau fc Mandidl2008l ; lMarassi et al.ll2009l) 
The detection of these astrophysical GWBs can provide in- 
sights into the cosmic star formation history and constrain 
some of the physical properties of compact objects, white 
dwarfs, neutron stars and black holes. Moreover, these sig- 
nals may act as foreground noise for the detection of cosmo- 
logical GWBs over much of the accessible frequency spec- 
trum. 

In this paper we consider the GWB produced by mag- 
netars, i.e. neutron stars with extremely high magnetic 
fields. Two classes of sources of high energy radiation in 
our Galaxy, the soft gamma repeaters (SGRs) and anoma- 
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lous X-ray pulsars (AXPs), are believed to host magnetars, 
that power their emission through the release of magnetic 
field energy. These two classes of high-energy sources share 
a number of features among which are the spin period in 
a fairly narrow and long range (in the 2 - 12 s range), the 
spin-down timescales (10 4 — 10 5 yr), the relatively faint per- 
sistent emission (typically 10 34 — 10 35 erg/s), the emission 
of sporadic short bursts (<C 1 s) wit h peak luminositie s in 
the 10 36 - 10 41 erg/s range (see e.g. iMereghettil 120081 and 
references therein). 

To successfully account for the observed features of 
both SGRs and AXPs, the magnetar model envisages that 
the neutron star possesses an internal magnetic field with 
strength B > 10 15 G which comprises both a toroidal and a 
poloidal component. The external B-field, on the contrary, 
is expected to be poloidal; its dipole strength is usually in- 
ferred to be in the 10 14 -10 15 G range, based on the observed 
spin-down rate (as well as other indications). Strong inter- 
nal magnetic fields (~ 10 15 to 10 16 G) will induce signifi- 
cant quadrupolar deformations in the neutron star structure; 
these may generate a detectable gravitational wave signal, 
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if t heir symmetry axis is not aligned with the spin axis (see 
e.g. ICutlerl ()2002h and references therein). 

Using population synthesis methods to evaluate the ini- 
tial period and the magnetic field di s tribut ions of magne- 
tars, iRegimbau fc de Freitas Pachecol (|2006t ) computed the 
GWB due to a magnetar population, assuming a purely- 
poloidal magnetic field configuration both inside of the star 
and in the magnetosphere. They found that the largest sig- 
nal is obtained for a type I superconductor neutron star 
model, and that the resulting closure energy density peaks 
at Qgw ~ 10~ 9 around 1.2 kHz; this is well below the sen- 
sitivity of the first generation of detectors, but it might be 
an interesting target for future detect ors, e.g. the Einstein 
Telesc opfl as discussed in detail in Re gimbau fc Mandid 
J2008h . 

In the present work, we reconsider the GWB gener- 
ated by magnetars, by using the cosmic star formation 
rate de nsity evolution predict ed by the numerical simula- 
tion of iTornatore et all l|2007h . and adopting several mag- 
netar models recently proposed in the literature. As a first 
example, we use the purely poloidal configur ations dis- 
cussed in IRegimbau fc de Freitas Pachecol (|2006l ) . We then 
consider tw i sted torus c onfigu rations, recent ly discussed in 
ICiolfi et al.l [|2009l . l2010h and in lFerraril (|20ld ): in these mod- 
els the poloidal component of the magnetic field extends 
throughout the star and in the magnetosphere, whereas the 
toroidal component is confined to a torus-shap ed region in- 
side t he star. Fi nally, we consider the model in lStella et al.l 
l|2005l) (see also iDall'Osso et all 120091 1. namely an internal 
field configuration dominated by the toroidal component 
with strength ~ 2 • 10 16 G (core-averaged value), and a 
poloidal field of ordinary strength (10 14 — 10 15 G). Our main 
purpose here is to assess how the uncertainties related to the 
internal magnetic field strength and its configuration affect 
the resulting GW signal and whether next-generation detec- 
tors will have the potential to reveal such a signal and shed 
light on the properties of magnetars. 

The plan of the paper is as follows. In Section [2] 
we briefly descr i be the numerical simulation performed by 
ITornatore et al. I (120071 ) that we use to predict the cosmic 
star formation rate evolution and the corresponding magne- 
tar birthrate. In Section[3]we sketch out the main features of 
the single source spectrum which describe the gravitational 
emission of a single magnetar, and introduce the different 
magnetar models considered in this work. In Section [3] we 
present the resulting density parameter of the GWB, flaw, 
and discuss its detectability. In Section [5] we discuss the ef- 
fect of the wobble angle on the generated background and its 
relevance for the detection of the signal. Finally, in Section 
[6] we draw our conclusions. 

Throughout our work we have adopt a ACDM cosmo- 
logical model with parameters Qm = 0.26, £7a = 0.74, 



THE MAGNETARS BIRTH RATE 
EVOLUTION 



WMAP results (|Spergel et al 



greemei 
Jill). 



Following M arassi et all i|2009h . we use the cosmic star for- 
mation rate den sity evolution predicte d by the numeri- 
cal simulation of ITornatore et all (|2007l ). For the present 
study, we consider the formation rate of Population II stars 
only; for these we adopt a Salpeter Initial Mass Function 
(IMF) $(M) oc M- (1+x) with x = 1.35 (normalized be- 
tween 0.1 — lOOM©), in regions of the Universe which have 
been already polluted by the first metals and dust grains 
ilSchneider et al.ll2002l 120031; lOmukai et al.ll2005h . 

Regimbau fc de Freitas Pachecol ( 200ll ) and 

iPopov et air T2009 ) derived the statistical properties 



of highly magnetized neutron stars (B^ 10 G) by using 
population synthesis methods; they showed that neutron 
stars born as magnetars represent 8-10% of the total 
simulated population of neutron stars. Here we assume a 
fraction fMNS = 10%; we further assume that magnetar pro- 
genitors have masses in the 8Mq — 40Mq range. It should 
be noted that the mass range of magnetar progenitors is 
still debated (see for in stance iFerrario fc Wickramasinghg 
2008; D avics et al.l f20091. However, the range we consider 
is sufficiently large to include the proposed evolutionary 
models. 

The top panel of Fig. [T] shows the redshift evolution 
of the cosmic star formation rate density inferred from the 
simulatiorQ. 

The number of magnetars formed per unit time out to a 
given redshift 2 can be computed by integrating the cosmic 
star formation rate density, p*(z), over the comoving volume 
element, while restricting the integral over the stellar IMF 
in the proper range of progenitor masses; that is 

AmnsW^/mns dz '-r, (, 2 \\ \ dM$(M), 
J az (i + z ) J 8M& 

(1) 

where the factor (1 + z) at the denominator accounts for the 
time-dilation effect, and the comoving volume element can 
be expressed as 

dV = 47rr 2 (j^J e(z)dz (2) 

e(z) = [n M (l + zf + fi A ]~* • 
The result is shown in the bottom panel of Fig. Q] 



3 MAGNETARS AS GRAVITATIONAL WAVE 
SOURCES 

The gravitational wave energy spectrum emitted by a single 
source is 

-1 



dEow 



= Eqw 



dt 



(3) 



where / e is the emission frequency. The gravitational wave 
luminosity of a neutron star with spin axis forming a wobble 



2 The results shown in Fig. [T] refer to the fiducial 



run m 



Tornatorc et al. I ||2007| 1 with a box of comoving size L 
1 http://www.et-gw.eu/ 10h~ 1 Mpc and N p = 2 X 256 3 (dark+baryonic) particles. 
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Figure 1. Top panel: redshift evolution of the comoving star 
formation rate density. Bottom panel: redshift evolution of the 
number of magnetars formed per unit time. 



angle a with the magnetic axis is composed of two contribu- 
tions, one at the spin frequency vr, one at its double 2vr\ 
it can be written as 



E G v 



2G 

5c 5 



T 2 2 6-2 / 2 , i a ■ 2 \ 

I € B uj sm a(cos a + losm a) 



(4) 



where e B is the star quadrupolar ellipticity induced by the 
magnetic field, cj is the angular velocity and I is the mo- 
ment of inertia. The term sin 2 a cos 2 a is from the emission 
component at the spin frequency, the term 16 sin 4 a from 
the component at twice the spin frequency. 

In strongly magnetized neutron stars the quadrupolar 
deformation is determined essentially by the magnetic field 
configuration and strength Moreover, it depends on the 
equation of state of matter in the stellar interior. 

The star loses rotational energy mainly due to electro- 
magnetic radiation and gravitational wave e mission (we shall 
negle ct other effects, e.g. relativistic winds; iDaH'Osso et all 
120091 ') . According to the well known vacuum dipole radiation 
model, the energy loss rate due to dipole radiation is given 



3 Fast rotation will also induce a non-negligible deformation, 
which, being symmetric with respect to the spin axis, does not 
contribute to the gravitational wave emission. 



\ E ROt\ = Q~^~ U S111 °' 



(5) 



where B p is the field strength at the magnetic poles, and R 
the stellar radius. The total spin-down rate obtained from 
Eqs. {4} and (0) is 



I • dip I . I . GW i 

1 BpR 6 a . 7 



(6) 



2G 



6 Jc 3 



cj" sm aH -/egoism a(l + 15sin a) 

5c 5 



Using the above quantity and remembering that the first 
term of the gravitational wave luminosity given in Eq. Q 
is emitted at f e = vr — u:/2n, while the second at f e = 
2ur, we can compute the single source emission spectrum 
according to Eq. Q. 

The choice of the initial spin period Po for the neu- 
tron star population sets an upper limit on the frequency 
ranges where the two components of the GW emission con- 
tribute to the gravitational wave background: the emission 
at vr contributes to frequencies below 1/Po, that at 2vr to 
frequencies below 2/Po- Therefore, for f e < 1/Po the gravi- 
tational wave background has both contributions, while for 
1/Po < f e < 2/Pq the only contribution comes from the 
emission at 2vr. Hence, the terms to be considered when 
computing the spectral energy density of the background 
are: 

for f e < ^ , 

dEcw 327T 4 G 2 2 x 3 
-I tBfe 



dfe 



5c 5 

BpR e 8n 2 G 2 ,2/, , , , . 2 \ 
— h „ 5 Ie B f e (l + 15 sm a) 



61c 3 



be 



+ sin q 



+ -Ie B fe(l + 15 sm a) 



6Jc 3 



for ^ < fe < £ , 
dEaw 32tt 4 G 



dfe 



for fe > % , 
dEcw 



5c 5 

BpR 6 
6/c 3 







5c 5 



T 2 2 j.3 



(7) 



. 27T G 2 f 2/, , -, c ■ 2 s 

+ 5c5 ie B / e (l + 15sm a) 



(8) 



If we assume B p = 10 14 - 10 15 G, R ~ 10 km, 7 ~ 10 45 
g cm 2 and f e ^ 1 kHz, we see that even for quadrupole ellip- 

B 2 R G 

ticities as large as 10 -4 the term -gj^r is much larger than 
8 5 C » , G /eg/ 2 ; in this case the contribution of gravitational 
wave emission to the spin-down is negligible. As shown in 
the following, this holds in most of the cases we consider. It 

B P R6 8tt 2 G r 2 f 2 



is worth noting that, when a^O and 
for / e < 1/Po the dominant term in Eq. © is 



dEGw 32-7r 4 G r2 2 j.3 

"J £fl/e 



5c r> 



P 2 p b 

6/c 3 



(9) 
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which does not depend on the wobble angle aQ. However, 



for -5- < / e < the dominant term is 



cLEgw 32n G 2 2 j.3 ■ 2 



67c 3 



(10) 



and it depends on a. 

It should be stressed that in general the wobble angle 
depends on time. The misalignment of magnetic and rota- 
tion axes causes, in the neutron star frame, the free preces- 
sion of the angular velocity around the magnetic axi s with 
period P vr<LC ~ P/|e_e|, where P is the spin period (|Jones! 
120021 ; IJones fc Anderssonll2002f ). The star internal viscosity 
damps such precessional motion and reduces the wobble an- 
gle towards the aligned configuration (a — 0), if the star has 
an oblate shape (es > 0), whereas it increases a towards the 
orthogonal configur ation a = n/2 ("spin-flip" ), if the shape 
is prolate (e B < 0) (I Jones! Il976l ; ICutled [20021 ). The second 
case is more favourable for gravitational wave emission. 

The timescale of the process is given by r a — nPo/tB, 
where Po/es is the initial precession period and n is the 
expected number of precession cycles in which the process 
takes place; estimates for slowly rotating neut ron stars indi- 
cate that n ~ 10 2 - 10 4 (|Alpar fe Sauls 1988), however the 
value of n is actually unknown (Cutler 2002). The evolution 
of the misalignment angle is relevant for our analysis only if 
the associated timescale, r a , is short compared to the spin- 
down timescale, r S£ j; conversely, if r a 2> r 8C i the process takes 
place when the source is no longer an efficient gravitational 
wave emitter. 

In principle, an accurate estimate of the GWB should 
account for (i) a proper distribution of the initial wobble 
angles for the magnetar population, and (ii) the evolution 
of such misalignment with time. This kind of analysis would, 
however, be affected by the wide uncertainties on both r Q 
and the initial angle distribution. Nevertheless, as shown 
in Section [SI we found that, from the point of view of the 
GWB detectability, the value of a and its eventual change 
with time do not significantly affect the results. Since our 
main interest is focused on detection prospects, we proceed 
here with the simplifying assumption that all magnetars are 
born with a = 7r/2 and that the misalignment evolution is 
ineffective. Then, in Section \S\ we will consider the effects 
of a generic wobble angle. 

An essential input for dEGw/df e is the magnetic field 
strength at the pole, B v : it determines the electromagnetic 
spin-down rate and, depending on the model, it may also 
affect the stellar deformations. To be representative of the 
entire population, the value of B p should be chosen as a 
suitable average. Such an average is uncertain at present; 
however the values of B p inf erred from AXPs and SGRs lie 
in the 10 14 - 10 15 G range (|Mereghettill200Sl ). Our choice 
here is to span this range by studying its two extremes, 



Br, 



10 14 and 10 la G. As we shall see, this translates into 



an uncertainty in the results that is negligible in comparison 
with the uncertainties associated to our poor knowledge of 
the internal field configuration. 



4 It should be no ted that if we remove th e term sin 2 a in Eq. (0 
(see for instance lOstriker fc Gunnl Il969f) , the a dependence in 



As we have seen, the overall GW emission de- 
pends on Pp, the initial spin period. Following 
iRegimbau fc de Freitas Pachecd (|200g) we set Po = 0.8 ms. 
For a = tt/2 this gives fT ax = 2/P = 2500 Hz. If a < tt/2, 
part of the gravitational wave emission is at the spin fre- 
quency and the corresponding contribution has a frequency 
cutoff f™ ax = I/-P0 = 1250 Hz. The chosen value of P 
implies a very fast spinning newborn neutron star, but still 
consistent with the believed range of neutron star spin rates 
at birth. We remark that, according to current scenarios 
of magnetar formation, strongly magnetized neutron stars 
are those that are born with period s of the order of ms, 
much faster than ordinary pulsars l|Duncan fc Thompson! 
1992). At the end of Section |4. II we will sketch the effect of 
assuming lower initial spin frequencies. 

In the following, we compute the gravitational wave 
background assuming different magnetic field models. 



3.1 Purely poloidal magnetic field 

The first two field configurations we consider describe a 
strongly magnetized neutron star endowed with a purely 
poloidal magnetic field, a nd have been used to evaluate 
the co rresponding GWB in lRegimbau fc de Freitas Pachecol 
|2006l ). 

It is well known that poloidal fields tend to make the 
star oblate (es > 0), while toroidal fields deform the star in 
a prolate shape (es < 0). Therefore, in the case we consider 
in this Section, ellipticity is always positive. Here we com- 
pute the gravitational wave emission according to Equations 
El) an d ((8|, setting, as in IRegimbau fc de Freitas Pachecol 
(2006), a — it/2 and assuming that the viscous evolution of 
the wobble angle (which tends towards a — in the case of 
positive ellipticities, thus reducing GW emission) is slower 
than the spin-down of the star. For a given poloidal configu- 
ration of the internal B-field, the GW output is maximized 
under the above assumption. 

The numerical inp uts to compute cLEg w I 'df e are B p , 
e_B, I and R. Following iKonno et al.l |2000), we write the 
quadrupolar ellipticity as 



£b = g 



73 2 P 4 
GhP 



(11) 



Eq. l|9jl is preserved. 



where M is the mass of the star and the value of 
the dimensionless (deformation) parameter g accounts 
for the magnetic field geomet r y and the EOS. As in 
IRegimbau fc de Freitas Pachecd (|2006l ). we consider two 
models with g = 13 (Model A) and g = 520 (Model B), re- 
spectively. The first model refers to an incomp ressible fluid 
star and a dipolar magnetic field l|Ferrardl 19541 ); similar val- 
ues are obtained in relativistic models based on polytropic 
equations of state l|Konno et al J 12000). Model B describes 
a scenario in which the neutron star core is a supercon- 
ductor of type I, implying that the internal magnetic field 
is con fined to the crustal layers (|Bonazzola fc GourgoulhorJ 
1 1996 ] ). This scenario give s much stronger deformations (see 
also fcolaiuda et al.ll200sl ). 

As previously discussed, we adopt two different values 
of B p : 10 14 and 10 15 G. The other parameters are fixed as 
follows: 7? = 10 km, M = 1.4 M Q and I = 10 45 g cm 2 . 
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3.2 Twisted torus configurations 

In order to account for the observed features of AXPs 
and SGRs, the magnetar model envisages that the internal 
magnetic field is a mixture of poloidal and toroidal com- 
ponents ( Duncan fe Thompson! Il992l ; iThompson fc Duncan! 
1 19931 . 1 19951 ). A strong internal toroidal components is ex- 
pected to form as a result of differential rotation shortly af- 
ter the birth of neutron star. The magnetic energy stored in 
this component provides the energy reservoir to power the 
bulk of the magnetar emission throu ghout its lifetime, in- 
cludi ng short bursts and giant flares ( Woods fe Thompson! 
2006). The poloidal component of the B-field (or at least 
part of it) extends outside the star, giving rise to a magneto- 
sphere. However, the detailed configuration of the magnetic 
field is presently unknown. Recent studies of the evolution 
of strongly magnetized stars in Newtonian gravity indicate 
that a particular magnetic field configuration, the so-called 
twisted torus, is a quite generic outcome of dynamical simu- 
lations and, due to magnetic helicity conservation, appears 
to be stable on dynamical time-scales jBraithwaite fc Spruit! 
12004 l200d : iBraithwaite fc Nordlundl l2006h . In this model, 
the poloidal field extends throughout the star and the exte- 
rior, whereas the toroidal field is confined in a torus-shaped 
region inside the star. 

Here we consider twisted torus equilibrium configura- 
tions of strongly ma gnetized neutron stars in the framework 
of general relativity (|Ciolfi et alj|2009l l20ldY The magnetic 
field includes contributions from higher multipoles (I > 1) 
coupled to the dipolar (I — 1) field which is usually assumed; 
an argument of minimal energy is adopted in order to estab- 
lish the relative weights of the different multipoles as well as 
the relative strength of toroidal and poloidal fields. In this 
model the poloidal fields dominate over the toroidal ones. 
As a consequence , these stars have always an oblate shape 
jCiolfi et al.ll2009l . l20ld ; lLander fc Jonesil2009l ). In order to 
account for the dependence on the equation of state, two dif- 
ferent EOSs are employed, named AP R2 and GNH3, wh ich 
span a realistic range of compactness (jCiolfi et al.ll20l5T ). 

As in the case of purely poloidal fields, we compute the 
gravitational wave emission spectrum according to Equa- 
tions (0, (JH]) and assuming a — tt/2. For a given magnetic 
field strength (B p = 10 14 and 10 15 G) and stellar mass 
(M = 1.4 Mq), the model provides the ellipticity, radius 
and moment of inertia of the neutron star; we have 



A: 



10 15 G 



■ 10" 



(12) 



with k = 9 (4), R = 14.19 (11.58) km, and / = 1.82 (1.33) 
10 45 gem 2 , for the EOS GNH3 (APR2). 



3.3 Toroidal-dominated magnetic field model 

The last model we consider is based on the hypotesis of very 
strong toroidal fields inside the star. We assume a magnetic 
field configuration with an internal toroidal field of ~ 2 • 10 16 
G (core-averaged value) , in addit ion to a polo i dal fie ld of or- 
dinary strength (10 14 - 10 15 Gl. lStella et~al1 (|2005l ) showed 
that toroidal field strengths of this order are needed to ex- 
plain the time-integrated emission of magnetars as inferred 
from the extremely bright giant flare that took place on 2004 
December 27 from SGR 1806-20 (which liberated an energy 



of 5 x 10 46 erg). Giant flares of this magnitude could result 
from large-scale rearrangements of t he core magnetic field 
or in s tabilities in the magnetosphere ([Thompson fc Duncan! 
l200ll ; iLvutikovl 120031 ). Such a huge toroidal field would in- 
duce prolate deformations as strong as e_e — —6.4 ■ 10~ 4 
([Cutler! I2002T ). Consequently, a newly-born fast spinning 
magnetar is expected to emit a strong gravitational signal 
whose frequency, initially in the 0.5 — 2 kHz range, decreases 
over a timescale of days (and whose strain correspondingly 
decreases too). This signal should be observable with the 
Advanced Virgo /L IGO class detectors up the distance of 
the Virgo cluster (|DaH'Osso et all 2009). The deformation 
associated to poloidal fields, whose strength is fixed by the 
choice of B p , tends to oppose the above deformation; how- 
ever in the cases considered here (B p = 10 14 and 10 15 G) 
the corresponding correction is negligible (of the order of 
10~ 4 — 10 -2 , respectively). The physical inputs we use in 
Equations (JTJl and (|SJ, in addition to a = 7r/2, en and B p , 
are R = 10 km and / = 10 45 gcm 2 . 

The gravitational wave emission predicted by the 
present model with B p — 10 14 G could be regarded as 
an upper limit among the different magnetar models (ex- 
cluding exotic scenarios), as in this case there occurs the 
most favorable combination of magnetic fields: an extremely 
strong toroidal field dominates the deformation while the 
lower value of the poloidal field strength results in a slower 
spin-down. In addition, the shape of the star is prolate and 
the evolution of the wobble angle a leads the axis of the 
magnetically-induced deformation towards the orthogonal 
configuration, resulting in stronger gravitational wave emis- 



3.4 Gravitational wave emission spectrum 

We present here the gravitational wave spectrum emitted by 
the models illustrated above. In Fig.[5]we plot cLEgw / df e as 
a function of the emitted frequency f e for the two purely 
poloidal models A and B (hereafter P-A and P-B), the 
twisted torus model (hereafter TT) with the two equa- 
tions of state considered (APR2, GNH3), and the toroidal- 
dominated model (hereafter TD). In the left (right) panel we 
assume B p — 10 14 G (B p = 10 15 G); as already discussed, 
these values define a likely range for B p . 

The first important indication which emerges from Fig. 
[5] is that the uncertainty related to the different magnetar 
models is always much higher (3-5 orders of magnitude) than 
the spread associated to the adopted range of B p . Let us now 
focus on the B p — 10 14 G case (left panel). The TD model 
is by far the most favorable for gravitational wave emis- 
sion, having the optimal combination of strong deformation 
and slow electromagnetic spin-down. The second strongest 
emission is obtained with the P-B model, where large defor- 
mations are achieved even for this lower field strength. The 
emission predicted by the P-A model is lower by more than 
three orders of magnitude, due to the difference in the g 2 fac- 
tor appearing in the expression of dEcw /df e - The two TT 
models are expected to give even weaker signals; they differ 
for the assumed EOS and, as expected, the one which gives 
less (more) compact stars, GNH3 (APR2), is associated to 
stronger (weaker) deformations and gravitational wave emis- 
sion. 

If we consider higher external poloidal fields (B p — 10 15 
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B p =10 14 G 



TT-APR2 — 

TT-GNH3 

P-A 

P-B 

TD 




f e [Hz] 
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Figure 2. Spectral gravitational energy emitted by a single source as a function of the emitted frequency, for the different models we 
consider: P-A and P-B stand for the purely poloidal models A and B, TT are the twisted torus model predictions for the two EOS 
considered (APR2 and GNH3), TD indicates the toroidal-dominated model. Left panel: B p = 10 14 G; right panel: B p = 10 15 G. 



G, right panel) the picture changes. For the P-A, P-B and 
TT models the value of B v controls both the gravitational 
wave luminosity, which scales as Bp (e_e oc Bp), and the 
spin-down rate, which has the electromagnetic contribution 
proportional to Bp 1 plus a very small correction due to grav- 
itational wave emission. As a result, dEcw I dfe oc Bp, which 
translates in a factor 100 increase from Bp = 10 14 to 10 15 G. 
Conversely, in the TD model the deformation is determined 
by the dominant toroidal field in the stellar interior (with 
poloidal field corrections up to ~ 1% for 10 15 G), and an 
increase in B p only results in a higher electromagnetic spin- 
down and in a smaller overall gravitational wave emission. 
As long as the electromagnetic spin-down dominates over 
the gravitational wave spin-down, dEGw/dfe is reduced by 
a factor of 100 from B p = 10 14 to 10 15 G. The final result 
is that when B p = 10 15 G, the prediction of TD and P-B 
models are comparable. 

It is worth noting that in all the considered models 
the contribution given by the gravitational wave emission 
to the spin-down is negligible, with the exception of the 
early time evolution in the TD model with B p = 10 14 G 
(hereafter TD14). This is shown in Fig. [2] where the en- 
ergy spectra are linear in logarithmic scale, reflecting the 
behaviour dEcw/df e oc /f, while the TD14 model is char- 
acterized by a lower emission level at high frequency, due to 
a non-negligible gravitational wave spin-down. This effect is 
even more evident in Fig. where we compare the grav- 
itational wave spectrum for the same TD model shown in 
the left panel of Fig. [2] with a TD model where gravitational 
wave spin-down is neglected (dashed line). It is clear that 
the gravitational wave contribution starts to be relevant at 
f e ~ 300 Hz. For all the other models we have discussed, 
this contribution becomes relevant at much higher emission 
frequencies. 



4 GRAVITATIONAL WAVE BACKGROUND 
FROM MAGNETARS 

In this Section we compute the GWBs produced by the dif- 
ferent magnetar models presented in the previous Section. 




f e [Hz] 

Figure 3. Spectral gravitational energy emitted by a single 
source according to the toroidal-dominated model as a function 
of the emitted frequency. TD-I is the same as TD in the left panel 
of Fig. [2] in TD-II (dashed line) the contribution of gravitational 
wave emission to the spin-down is neglected. 



Following Mar assi et al.l l|2009t ). the spectral energy density 
of the GWB can be written as 

dR(M,z){-^}, (13) 



dE 



dSdfdt 

where dR(M, z) is the differential source formation rate 
p*(z) dV 



dR(M, z) 



(1 + z) dz 



$(M)dMdz 



(14) 



and ( ) is the locally measured average energy flux emit- 
ted by a source at distance r. For sources at redshift 2 it 
becomes 

{l + zf dE G w, 



dE 



-[/(! + *) 



(15) 



-dSdf 4nd L (z) 2 dfe 

where / = / e (l + z) - is the redshifted emission frequency 
f e , and dh{z) is the luminosity distance to the source. 

It is customary to describe the GWB by a dimen- 
sionless quantity, the closure energy density Qgw(/) = 
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Figure 4. The predicted closure energy density (J2qw) ^ a function of the observational frequency, for the different magnetar models 
discussed in this paper: P-A and P-B stand for purely poloidal models A and B, TT are the twisted torus model predictions for the 
two EOS considered (APR2 and GNH3), TD is the toroidal-dominated model. Left panel: B p = 10 14 G; right panel: B p = 10 15 G. In 
both panels the shaded region indicates the foreseen sensitivity of the Einstein Telescope, and the horizontal dotted line (CGWB) is 
the upper limit on primordial backgrounds generated during the Inflationary epoch. A given background is detectable by the Einstein 
Telescope if the corresponding signal-to-noise ratio is larger than the detection threshold S/N = 2.56 (see text). 



p cr 1 (dp gw /dlogf), which is related to the spectral energy 
density by the equation 



!iow(/) = 



(IE 



dSdfdt 



(16) 



where p cr = 3Hq/8itG is the cosmic critical density. 

In Fig. [H we show Qgw as a function of the observa- 
tional frequency for the different magnetar models. Differ- 
ences in the predicted stochastic backgrounds reflect differ- 
ences in the corresponding single magnetar emission spec- 
trum. The maximum amplitude is always achieved around 1 
kHz: in the left panel, it ranges from ~ 4- ltr 16 to ~ 2- 10 -8 , 
while in the right panel the range is ~ 4-10~ 14 -2-10~ 9 . The 
higher value is obtained with the TD model in the first case, 
and with the P-B model in the second case; the lower value is 
given in both cases by the TT-APR2 model. In both panels, 
model predictions are compared with the foreseen sensitiv- 
ity of the Einstein Telescope (shaded region) and with the 
upper limit to primordial backgrounds generated during the 
Inflationary epoch (horizontal dotted line labelled C GWB). 
The la tter contribution is estimated from Eq. (6) of I Turner] 
i|l997l ) assuming a tensor/scalar ratio of r = 0. 3 and no run- 
ning spectr al index of tensor perturbations (jKinnev et al.l 

l2ood . l2ooaf) . 

It is clear from the figure that the background generated 
by magnetars in the kHz region is higher than the primordial 
background, independently of the specific magnetar model 
considered; in addition, specific magnetar models lead to 
a cumulative signal which is potentially detectable by the 
Einstein Telescope. A more quantitative assessment of the 
detectability is reported in the following Section. 



4.1 Detectability 

The gravitational wave signal produced by the magnetar 
population can be treated as continuous. Indeed, if Ar gw is 
the average time duration of a signal produced by a single 
magnetar, and dR(z) is the number of sources formed per 



unit time at redshift z, the duty cycle D out to redshift z, 
defined as 



D(z) 



■ f dR(z)Ar gw (l + z) 
Jo 



(17) 



satisfies the condition^ D>1. Consequently, the stochastic 
signal appears in the detector outputs as a time-series noise 
which, by the central limit theorem, is expected to have 
a Gaussian-normal distribution function. In this case, as 
sugges ted by I Allen fc Romanol (|l999T ): [Regimbau fc Mandicl 
(2008), the optimal detection strategy is to cross-correlate 
the output of two (or more) detectors, assumed to have in- 
dependent spectral noises. 

The optimized S/N fo r an integration time T is given 
by I Allen fc Romanol l(l999l ). 



9#o 4 
50tt 4 



df 



7 2 (/)"gw(/) 
/ 6 A(/)P 2 (/) ' 



(18) 



where Pi(f) and Pz(f) are the power spectral noise densities 
of the two detectors, and 7 is the normalized overlap reduc- 
tion function, characterizing the loss of sensitivity due to 
the separation and the relative orientation of the detectors. 

The sensitivity of detector pairs is given in terms of the 
minimum detectable amplitude for a flat spectrum J7min 
(fiiviiN = const) defined as 



1 IOtt 
•(erfc~ 1 (2a 



df- 



7 2 (/) 



/ 6 A(/)P 2 (/) 
-erfc- 1 (2 7 )), 



-1/2 



(19) 



where T is the observation time, a the false alarm rate, 7 the 
detection rate and e rfc~ the complementar y error function 
(for more details see lAllen fc Roniandll999M . 

If we consider the cross-correlation of two detectors with 



5 If we take At 



T sd we have always D higher than 10 
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Figure 5. The TD14 background is plotted for different wobble 
angles, spanning the range 7r/60 — vr/2, and compared with the 
Einstein Telescope sensitivity. 



the sensitivity of the Einstein Telescope (T. Regimbau, pri- 
vate communication), we get Omin = 1.13 • 10 -11 for an in- 
tegration time T of one year, a false alarm rate a = 10% and 
a detection rate 7 = 90%; these values, inserted in Eq. (|18[) . 
lead to a detection threshold S/N of 2.56. 

A given background is detectable by the Einstein Tele- 
scope if the corresponding signal-to-noise ratio given by 
Eq. (|18|) is larger than the detection threshold. For instance, 
the predicted f^GW for the purely poloidal model P-B with 
B p = 10 15 G gives S/N = 2.49, that is slightly smaller 
than such threshold; consequently, there is no chance to de- 
tect this signal. Conversely, in the most optimistic magnetar 
model TD14 (toroidal-dominated model with B p = 10 14 G) 
we obtain S/N — 40, a very promising value. This result 
leads to the conclusion that third-generation gravitational 
wave detectors, such as the Einstein Telescope, hold the po- 
tential to reveal the cumulative GW signal from magnetars 
in the universe. 

It is worth noting that the above results refer to the 
assumed initial spin period of Po = 0.8 ms. A higher value 
would lead to a lower frequency cutoff and, consequently, 
to a weaker GWB. For the TDn model, for example, the 
detection threshold S/N — 2.56 corresponds to Po = 5.2 
ms. Hence the GWB would still be detectable up to this 
value. 



5 WOBBLE ANGLE EFFECTS 

So far we have assumed a constant misalignment a = 7t/2 
between the spin and the magnetic axis, in which case the 
GW signal is emitted only at twice the spin frequency f e = 
2pr = uj/n. For a generic misalignment, we have also the 
emission at the spin frequency. 

We now focus on the model TD14 and explore the con- 
sequences of q < 7r/2. In this model the stellar deformation 
induced by the magnetic field are larger; being the most 
optimistic model for gravitational wave emission, this case 
allows to clearly show the effects of the wobble angle on 
detectability. 

In Fig. [5] we compare the sensitivity of the Einstein 
Telescope with the background generated by model TD14 



for different (constant) values of the wobble angle. The plot 
clearly shows that when a > 7t/4 there is a single dom- 
inant contribution with the frequency cutoff at 2500 Hz, 
while for smaller angles there is a dominant contribution 
with f max — 1250 Hz and a secondary contribution with 
lower amplitude extending up to f max — 2500 Hz. Similar 
effects hold for the alternative magnetar models which have 
been presented in the previous Sections. A difference be- 
tween the models potentially arises if the timescale for the 
evolution of the wobble angle is short compared to the spin- 
down timescale (see Section [3]): in this case, the star rapidly 
tends (i) to the orthogonal configuration for the TD model, 
thus increasing gravitational wave emission, and (ii) to the 
aligned configuration, for models P-A, P-B and TT models, 
thus decreasing the emission. 

As shown in Fig. \S\ for the TD14 model (as well as for 
the other models) the gravitational wave backgrounds corre- 
sponding to different wobble angles exhibit significant differ- 
ences at large frequencies, approximately above ~ 800 Hz, 
where the Einstein Telescope sensitivity is too low even 
for this model; therefore, the signal detectability is only 
marginally affected. Variations in the signal-to-noise ratio 
are at most 2-3% in the TD14 case, and if the gravitational 
wave background is weaker (e.g. for B p > 10 14 G) the effects 
on the S/N are even smaller. 

We can conclude that the initial value of a and its evo- 
lution in time do not have significant effects on the GWB 
detectability with the Einstein Telescope. 



6 CONCLUSIONS 

In this paper, we estimated the GWB produced by magne- 
tars. We used a cosmic star formatio n history obtained from 
a numerical simulation performed bv lTornatore et all (|2007T l 
and assume that 10% of stellar progenitors with masses in 
the range 8 — 40M© lead to magnetars with magnetic field 
strength in the 10 14 - 10 15 G range. 

Since our present understanding of the physical prop- 
erties of magnetars is still poor (the internal structure of 
the magnetic field and initial spin frequency being among 
the major uncertainties), we have explored the consequences 
for gravitational wave emission of different magnetar models 
proposed in the literature. 

Our analyis shows that different models produce a 
spread in the resulting gravitational wave emission which 
is much higher than that produced by adopting different 
values for the magnetic field strength. In particular, we find 
that: 

• Toroidal-dominated models, with an internal toroidal 
field of ~ 2 • 10 16 G and an external poloidal field of 10 14 G, 
proposed by Stella et al. (2005) to explain the 2004 giant 
flare from SGR 1806-20, generate the largest gravitational 
wave background, which could be detected in the frequency 
range between ~ 50 and ~ 600 Hz by third generation grav- 
itational wave detectors such as the Einstein Telescope. Us- 
ing correlated analysis of Einstein Telescope outputs, the 
estimated signal-to-noise ratios could be as high as 40. 

• When larger poloidal fields, 10 15 G, are considered, the 
largest gravitational background is generated by magnetar 
models with purely poloidal fields, and a superconductor 
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type-I core; in this case, the internal magnetic field is con- 
fined to the crustal layers, leading to strong deformations. 
Since deformations are produced by the internal toroidal 
field, Toroidal-dominated models are less effective because 
the increase in the poloidal field strength leads to a higher 
electromagnetic spin-down and to a lower gravitational wave 
emission. 

• A comparison between the estimated magnetar GWB 
and the upper limit to the primordial background predicted 
by Inflationary scenarios (the horizontal dotted line in Fig. 
[4] labelled CGWB) shows that, for the models of magnetar 
we consider, the magnetar GWB is always larger than the 
primordial background in some region of frequency (the only 
exception is the twisted-torus model TT-APR2 with B p = 
10 15 G). For instance, for the toroidal dominated models 
TD15 and TD14 this is true, respectively, for / > 13 Hz and 
/ > 4 Hz. Thus, the GWB generated by magnetars may 
act as a limiting foreground for the future detection of the 
primordial background even at frequencies as low as few tens 
of Hz. 

We have also investigated the consequences on the re- 
sulting gravitational wave background of different values for 
the misalignement angle between the rotation and magnetic 
field axes. We find the largest effects to be at high frequen- 
cies, above ~ 800 Hz; thus the detectability of the largest 
backgrounds with the Einstein Telescope is only marginally 
affected, with fractional variations of the signal-to-noise ra- 
tios of at most 2-3%. 
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